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In the past few decades, traditional integral transform methods such as Fourier and Laplace trans¬ 
forms have commonly been used to solve engineering problems. These methods transform differential 
equations into algebraic equations which are easier to deal with. However, these integral transform 
methods are more complex and difficult when applying to nonlinear problems. 

The HPM, proposed first by He[2,3], for solving the differential and integral equations, linear and 
nonlinear, has been the subject of extensive analytical and numerical studies. The method, which is 
a coupling of the traditional perturbation method and homotopy in topology, deforms continuously 
to a simple problem which is easily solved. This method, which does not require a small parameter 
in an equation, has a significant advantage in that it provides an analytical approximate solution to 
a wide range of nonlinear problems in applied sciences. 

The HPM is applied to Volterra’s integro-differential equation [4], to nonlinear oscillators [5], bi¬ 
furcation of nonlinear problems [6], bifurcation of delay-differential equations [7], nonlinear wave 
equations [8], boundary value problems [9], quadratic Riccati differential equation of fractional or¬ 
der [1], and to other fields [10-18]. This HPM yields a very rapid convergence of the solution series 
in most cases, usually only a few iterations leading to very accurate solutions. Thus He’s HPM is a 
universal one which can solve various kinds of nonlinear equations. 

Adomian decomposition method(ADM), which was introduced by Adomian[19], is a semi numeri- 
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cal technique for solving linear and nonlinear differential equations by generating a funtional series 
solution in a very efficient manner. The method has many advantages: it solves the problem di¬ 
rectly without the need for linearization, perturbation, or any other transformation; it converges 
very rapidly and is highly accurate. 

Differential transform method (DTM),which was first applied in the engineering field by Zhou[8],has 
many advantages :it solves the problem directly without the need for linearization,pertuburation,or 
any other transformation.DTM is based on the Taylors series expansion.lt is different from the tradi¬ 
tional high order Taylor series method,which needs symbolic computation of the necessary derivatives 
of the data functions.Taylor series method computationally takes a long time for larger orders while 
DTM method reduces the size of the computational domain,without massive computations and re¬ 
strictive assumpions,and is easily applicable to various physical problems .The method and related 
theorems are well adressed in [9,10].Let us consider the (2+l)-dimensional Burgers’ equation with 
damping term as, 


under the initial condition 


ILt — U x Uy T 'll T Uxyi 


( 1 ) 


u(x,y,0) = u 0 (x,y). (2) 

Burger’s equation is the simplest second order NLPDE which balances the effect of nonlinear 
convection and the linear diffusion.In this work ,we have employed the Homotopy Perturbation 
Method(HPM),DTM and ADM to solve the (2+l)dimensional Burger’s equation with damping 
term. 


1.1 Homotopy Perturbation Method (HPM) 

To describe the HPM, consider the following general nonlinear differential equation 


A(u) - f(r) = 0,reD, 


(3) 


under the boundary condition 


-B(%|j^|= o ,re0f2, (4) 

Where A is a general differential operator, B is a boundary operator,f(r)is a known analytic function 
and dfl is a boundary of the domain fl.The operator A can be divided into two parts L and N, 
Where L is a linear opertor while N is a nonlinear operator.Then Equation(3)can be rewritten as 


L{u) + N(u) - /(r) = 0, (5) 

Using the homotopy tecnique, we construct a homotopy: 

V(r,p) : DX[0,1] R 
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Which satisfies 

H(V,p) = (1 -p)[L(V) - L(«o)] +P[A(V) - f(r)] 


or 


H(V,p) = L(V) - L(u 0 ) + pL(u 0 ) + p[N(V) - f(r)],reQ,pe[ 0,1], ( 6 ) 

Where pe[ 0,1] is an embedding parameter, uq is the initial approximation of Equation(3) which 
satisfies the boundary conditions.Obviously,considering Equation(6),we will have 

H(V,0) = L(V)-L(u 0 ) = 0, 


H(V, 1) = A(V) - f(r) = 0, 


(7) 


changing the process of p from zero to unity is just that V(r,p)from Uo(r) to u(r).In topology,this is 
called the deformation also A(V)-f(r) and L(u)are called as homotopy .The homotopy perturbation 
method uses the homotopy parameter p as an expanding parameter [23-25] to obtain 

OO 

V = v 0 +pv i +p 2 v 2 +p 3 v 3 + ... = ^2p n v n . (8) 

n—0 

p —>1 results the appropriate solution of equation(3)as 


u = lim V = Vo + i’i + v 2 + ■■ 
p-* i 


n=0 


(9) 


A comparision of like powers of p gives the solutions of various orders series(9)is convergent for most 
of the cases.However convergence rate depends on the nonlinear, N(V) He[25]suggested the following 
opinions: 

1. The second derivative of N(V)with respect to v must be small as the parameter p may be relatively 
large. 

dN 

2. The norm of ZA 1 —— must be smaller than one so that the series converges. In this section, we 

ou 

describe the above method by the following example to validate the efficiency of the HPM. 


Example: 1 

Consider the (2+l)-dimensional Burgers’ equation with darning as, 


under the initial condition 


U t — Vj X y T Ux'Uy T u, 


u{x, y, 0) = u 0 {x, y) = x + y. 

Applying the homotopy perturbation method to Equation(10),we have 


( 10 ) 


( 11 ) 



U t +P[(-U xy 

— uu x — zt)] = 0 


(12) 
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In the view of HPM, we use the homotopy parameter p to expand the solution 


u(x,y,t) = u 0 +pui + p 2 u 2 + ■■■ 


(13) 


The approximate solution can be obtained by taking p=l in equation (13)as 


u(x,y,t) = u 0 + ui + u 2 + ... (14) 

Now substituing equation(12)into equation(ll)and equating the terms with identical powers of p, 
we obtain the series of linear equations, which can be easily solved.First few linear equations are 
given as 


P 


o . duo 
' dt 


dui 


dt 


d 2 u 0 

dxdy 


du 0 du 0 , 

+ -w- q -h Uo ■ 

dx dy 


du 2 


dt 


d 2 Ul 

dxdy 


dui duo 
dx dy 


duo du i 
dx dy 


+ U\. 


Using the initial condition(ll), the soution of Equation(15)is given by 


u{x, y, 0) = u 0 (x, y) = (x + y). 


Then the solution of Equation(16)will be 


ui(x,y,t) 



duo duo 
dx dy 



dt. 


ui(x,y,t) = (x + y)t + t. 

Also, We can find the solution of Equation(17)by using the following formula 


t 

ui(x,y,t) = j 


( d 2 ^ 1 , duo duo dui . , 
\dxdy dx dy dx dy Ul 


f 2 q+2 

u 2 (x,y,t) = (x + y )— + — 

etc.therefore, from equation(18),the approximate solution of equation(10)is given as 

t 2 3t 2 

u{x, y,t) = (x + y) + {x + y)t + t+ — + — + ... 

Hence the exact solution can be expressed as 


u(x,y,t) = {x + y)e t + (n\ + l)(e‘ - 1), 


(15) 

(16) 

(17) 

(18) 

(19) 

( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 


provided that 0 < t < 1 
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1.2 Adomain Decomposition Method(ADM) 

Consider the following linear operatpor and their inverse operators: 

,-A, » 

*' 9t' dxdy' 

t xy 

Li 1 = J (• )dt,L x>y = j J(jdrd'y. 

0 0 0 

Using the above notations,Equation(l)becomes 

Lti^) = L X y{lL) U X Uy , 

Operating the inverse operators L^ 1 to equation(25)and using the initial condition gives 

u(x, y , t) = u 0 (x, y , t) + L t _1 (T Xi . y (u)) + Li 1 (u x u y ), (26) 

The decomposition method consists of representing the solution u(x,y,t)by the decomposition series 

OO 

u(x, y,t) = ^2 u q ( x > 0- (27) 

9 =o 

The nonlinear term u x u y is represented by a series of the so called Adomain polynomials, given by 

OO 

u = A q (x,y,t). (28) 

9=0 

The component u g (x,y,t ) of the solution u(x,y,t)is determined in a recursive manner.Replacing the 
decomposition series (27)and(26)gives 

OO OO 

^ 2 u g(x,y,t ) = u 0 (x,y,t) + Li 1 (L x , yt t(u)) + Li l ^A q (x,y,t) (29) 

9=0 9=0 

According to ADM the zero-th component uo(x,y,t) is identified from the initial or boundary cndi- 
tions and from the source terms.The remaining components of u[x , y , t ) are determined in a recursion 
manner as follows 


u 0 (x,y 7 t) = u 0 (x,y), 

(30) 

u k (x,y,t) = T t _1 (Tx, y (u)) + Li 1 (A k ),k > 0, 

(31) 


Where thhe adomain polynomials for the nonlinear term u x u y are derived from the following recursive 
formulation 


Ak 


1 d k 




,k 


A=0 


0 , 1 , 2 ,... 


First few Adomian polynomials are given as 


duo dui duo 

A) = u 0 -—, Ai = u 0 - -hui-—, 

OX i OX i OX i 


(32) 
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A du 0 
A 2 = u 2 w- 

OX i 


Ml 


dui 

dx\ 


u 0 


du 2 
dx\' 


using equation(31)for the adomain polynomials Ak,~we get 


(33) 


uo{x,y,t) = u 0 (x, y), (34) 

ui(x,y,t) = Ll l (L x ,y(u 0 )) +Lt 1 (A 0 ), (35) 

u 2 (x, y, t) = Li\L x , y { Ul )) + L t _1 (Ai), (36) 

and so on.Then the q-th term, u q can be determined from 

9-1 

u q = '^2 u k(x,y,t). (37) 

o 


Knowing the components of u, the analytical solution follows immediately. 


1.2.1 Computational illustrations of ADM for (n+l)-dimensional Burgers’ equation 
with damping 

using Equations(32) and (33), first few components of the decomposition series are given by 


uo(x,y,t) = (x + y ), 

(38) 

ui (x,y,t) = (x + y)t + t, 

(39) 

t 2 3 1 2 

u 2 (x,y,t) = (x + y)— + —, 

(40) 

t 3 7 1 3 

u 3 (x,y,t) = (a; -h 9 /) — + —, 

(41) 

OO 

u(x,y,t) = ^2Uk(x,y,t ) 
k =0 


■u 0 (x,y,t) + m (x,y,t) +u 2 {x,y,t) + ..., 



oj-2 42 43 743 

= {x + y) + {x + y)t + t + — + (x + y)— + (m + 9/) — + ~^~ + - ( 42 ) 

Hence, the exact solution can be expressed as 

u(x, y, t) = (x + y)e t + (n\ + l)(e* - 1) (43) 

1.3 Differential Transform Method(DTM) 

In this section,we give some basic definitions of the differential transformation.Let D denotes the 
differential transform operator and D -1 the inverse differential transform operator. 


@ IJTSRD 

Available Online @ www.ijtsrd.com 

Volume - 2 

Issue - 4 

May-Jun 2018 • Page: 2119 





International Journal of Trend in Scientific Research and Development (IJTSRD) ISSN: 2456-6470 


1.3.1 Basic Definition of DTM 
Definition:! 


If u{x 1 , 3 : 2 , ■■■,x n ,t) is analytic in the domain fl then its (n+1)- dimensional differential transform 
is given by 


U (fci . k 2 ,/c n , 


1 


fci!. fc 2 !, ■. ■, k n \, k n -\. i! 

-\-k2-\-...-\-kn-\-k n +i 


(44) 


r) kl r) k2 r) kn r)‘ 

DXi Ux2 • " U X n 67 


,k n +1 
x t 


.u(xi,x 2 , ...,x n ,t)\xi = 0 ,a ;2 = 0, ...,x n = 0,t = 0 


where 


OO OO OO OO 

u(x 1 ,X2,...,x n ,t) = U (ki,k 2 ,...,k n ,k n+1 ) .x^ 1 x l 2 2 ...x kn t kn + 1 

kl=0 &2 = 0 kn= 0 &n +1=0 (45) 

— -D [C/ (^1? ^ 25 •••> &n+l)] 


Definition^ 


Ifu(aji,x 2 , —,x n ,t) = D 1 [u(k 1 ,k 2 ,...,k n ,k n+ i)\,v(xi,x 2 ,...,x n ,t) = D 1 [V(k 1 , k 2 , k n , fc„+i)],and 
® denotes convolution,then the fundamental operations of the differential transform are expressed 
as follows: 


(a).D[u(x i, x 2 , x„, t)v(x i, x 2 , x n , f)] = U(ki,k 2 , ..., k„, k n+ 1) <g> V(k i, fc 2 , ■■■, k n , k n+1 ) 

fcl &2 fcn fcn+1 

= U(ai, k 2 — a 2 , ..., fc„+i — a n +i)V(ki — ai, a 2 ) a„+i) 


ai=0a2=0 a n =0 a n +i=0 


(46) 


(6).D[au(a;i,x 2 , ...,x n ,t) ± / 3 v(xi,x 2 , ...,x n ,t)] = aU(ki, J^2 5 •••5 &n+l) 

£)ki + Al2+. • .+fcn + &n-|-l 

(d J t t 2 » fc „ +1 ■»( a; i^2,-,a: n ,f) = (fcl + l)(fci +2)...(fci + n)(fc 2 + l)(fc 2 + 2) 

(7^2 •••C'x n Oxt 

...(fc 2 + r 2 )...(fc n+ i + l)(fc n _|_i + 2)...(fc n+1 + r n+ i).t/(fci +ri ,..., k n+ 1 + r„ + i) 


(47) 

(48) 


1.3.2 Computational illustrations of (2+l)-dimensional Burgers’ equation with damping 

Here we describe the method explained in the previous section, by the following to validate the 
efficiency of the DTM Consider the (2+l)-dimensional Burgers’ equation with damping as, 

u t = u xy + u x u y + u, (49) 

Subject to the initial condition 


u{x, y, 0) = u 0 (x, y) = x + y 


(50) 


Talking the differential transform of equation(49).,we have 

{ks + l)D(fci, fc 2 , fc 3 + 1) = (fci + 2)(fci + 1)C/(fci + l,fc 2 + 1, fc3) + f/(fci,fc 2 ,fc 3 ) 

+U (oi, fc 2 — a 2 , fc 3 — a 3 ) + U (fci — ai, a 2 , fc 3 — a 3 ) 


(51) 
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From the initial condition equation(50),it can beseen that 


OO OO 


i{x,y, 0) = U(k 1 ,k 2 ,0)(x kl .y k2 ) =x + y 

k±=0 /c2—0 


Where 


U(hi, k 2 ,0) = 


1 if ki = 1 , kj = 0 , i ^ j,i,j=l ,2 


0 otherwise 

Using eq.(53)into eq.(52)one can obtain the values of U(k\, k 2 , k 3 , £q) as follows 


U {k \, k 2 ,1) = 


1 if ki = 1 , kj = 0, i ^ j,i,j=l ,2 


0 otherwise 


U(ki,k 2 , 2) = 


1 if ki = 1 , kj = 0 j, i , j = 1,2 


0 otherwise 


Then from eqn.(45)we have 


OO OO 

u(x,y,t)= ^ u ( k h k 2 ,k 3 )x kl y k2 t k3 

k\ —0 /c2—0 


3 1 2 


. t 3 7 1 3 


= {x + y) + {x + y)t +1 + — + (x + y)— + (x + y)- + — 


Hence,the exact solution is given by 


(52) 


(53) 


(54) 


(55) 


(56) 


u(x , y, t) = (x + y)e + (n! + l)(e - 1 ) 


(57) 


1.4 Conclusion 

l.In this work,homotopy perturbation methhod,DTM and adomian decomposition method have 
been successfully applied for solving (2+l)-dimensional Burgers’ equation with damping term. 

2 . The solutions obtained by these methods are an infinite power series for an appropriate initial 
condition, which can,in turn be expressed in a closed form, the exact solution. 

3. The results reveal that the methods are very effective, convenient and quite accurate mathe¬ 
matical tools for solving the (2+l)-dimensional Burgers’ equation equation with damping. 

4. The solution is calculated in the form of the convergent power series with easily computable 
components. 

5. These method, which can be used without any need to complex computations. 
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